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ABSTRACT 

The high eccentricities of the known extrasolar planets remain largely unex- 
plained. We explore the possibility that eccentricities are excited in the outer 
parts of an extended planetary disk by encounters with stars passing at a few 
hundreds of AU. After the encounter, eccentricity disturbances propagate inward 
due to secular interactions in the disks, eventually exciting the innermost planets. 
We study how the inward propagation of eccentricity in planetary disks depends 
on the number and masses of the planets and spacing between them and on the 
overall surface-density distribution in the disk. The main governing factors are 
the large-scale surface-density distribution and the total size of the system. If the 
smeared-out surface density is approximated by a power- law E(r) oc r~ g , then ec- 
centricity disturbances propagate inward efficiently for flat density distributions 
with q < 1. If this condition is satisfied and the size of the planetary system is 
50 AU or larger, the typical eccentricities excited by this mechanism by field star 
encounters in the solar neighborhood over 5 Gyr are in the range 0.01-0.1. Higher 
eccentricities (>0.1) may be excited in planetary systems around stars that are 
formed in relatively dense, long-lived open clusters. Therefore, this mechanism 
may provide a natural way to excite the eccentricities of extrasolar planets. 

Subject headings: stars: planetary systems - planets and satellites: formation 



1. Introduction 

One of the remarkable features of the ~ 120 known extrasolar planets is their rela- 
tively high eccentricities 1 , most far larger than seen in the giant planets of the solar system. 
Many mechanisms to produce these high-eccentricity orbits have been offered. The most 
popular involve planet-planet gravitational scattering (Rasio & Ford 1996; Weidenschilling 



1 Princeton University Observatory, Princeton, New Jersey 08544 

1 A catalog of extrasolar planets is maintained by Paris Observatory (http://www.obspm.fr/planets) 
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k Marzani 1996; Lin k Ida 1997; Ford et al. 2001; Yu k Tremaine 2001; Chiang, Fischer, 
k Thommes 2002; Terquem k Papaloizou 2002), either through close encounters, physical 
collisions, secular interactions or resonant interactions. Other eccentricity excitation mech- 
anisms include interactions with the gaseous protoplanetary disk (Goldreich k Tremaine 
1980; Papaloizou et al. 2001; Goldreich k Sari 2003; Ogilvie k Lubow 2003), or interactions 
with the planetesimal disk (Murray et al. 1998). 

All of these proposed mechanisms have shortcomings. Planet-planet scattering through 
close encounters appears to create the wrong eccentricity distribution — there may be too 
many planets on low-eccentricity orbits, formed by physical collisions or tidal captures dur- 
ing the scattering process (Ford et al. 2001; Goldreich k Sari 2003, but see Ford, Rasio, k Yu 
2003). Discussions of eccentricity excitation through mean-motion resonances assume that 
planet-disk interactions induce evolution of semimajor axes leading to resonant capture, but 
neglect the effects of these interactions on eccentricity. Interactions with the gaseous plan- 
etesimal disk can either damp or excite eccentricity, depending sensitively on the formation 
of gaps at Lindblad resonances, nonlinear saturation of corotation resonances (Ogilvie k 
Lubow 2003), and the rate of viscous diffusion across resonances (Goldreich k Sari 2003), all 
of which are difficult to calculate reliably. Excitation through interactions with the planetes- 
imal disk may require more massive disks than are indicated by observations and may only 
be efficient for companions with masses much higher than those of most extrasolar planets 
(Papaloizou et al. 2001). 

There is strong and growing evidence that most forming stars are surrounded by ex- 
tended disks, with radii of hundreds of AU. This evidence includes CO observations of mas- 
sive gas disks around young stars, disks detected in Orion and other star-forming regions 
by optical emission from their photoionized surfaces, the infrared excess at wavelengths of 
25-100/i discovered by IRAS around bright nearby stars such as Vega, dust disks found op- 
tically around nearby stars such as (5 Pictoris, and imaging by the NICMOS camera on the 
Hubble Space Telescope of near-infrared radiation scattered by dusty disks around stars such 
as HR 4796 (see Koerner 2001 and Zuckerman 2001 for reviews). It may well be that the 
dust around older stars such as Vega and f3 Pictoris is "second-generation" dust, released in 
collisions between larger objects (planetesimals or comets) that were formed by the accretion 
of primordial or first-generation dust particles. Thus the absence of dust around older stars 
may arise either because the dust has been dispersed or destroyed or because it has been 
incorporated into larger bodies such as planets. 

Indirect evidence for planets at large distances from the host stars comes from the 
morphologies of the dust debris disks of Vega and e Eri. Disk asymmetries and dust concen- 
trations are often interpreted as resulting from the dynamical influence of an unseen massive 
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planet. Modeling suggests that the planetary companion has a semimajor axis a > 30 AU 
in the Vega system (Wilner et al. 2002) and about 55-65 AU in the e Eri system (Ozernoy 
et al. 2000). Formation of giant planets by core accretion is difficult at these distances, and 
it has been suggested that planets can migrate outward to distances of up to tens of AU 
through interactions with the gas disk in young systems (Veras & Armitage 2004). 

We propose to explore the possibility that the eccentricities of the extrasolar planets 
arise because of secular interactions in a long-lived extended planetary or protoplanetary 
disk. We suggest that eccentricities are excited in the outer part of the disk by a passing star 
and propagate inward through the disk, somewhat like a wave. We explore how the inward 
propagation of eccentricity in disks depends on the overall surface-density distribution. 

In Section 2 we review the definitions and the approach of secular perturbation theory. 
In Section 3 we describe the initial excitation of the eccentricity and in Section 4 we study 
the propagation of eccentricities via secular interactions. We discuss typical eccentricities 
produced by this mechanism in Section 5 and conclude in Section 6. 



2. Secular evolution of a planetary system 

In this Section we consider a planetary system around a star with mass M c . There are N 
planets with masses rrii <C M c and orbital elements a« (semimajor axes), Si (eccentricities), 
Wi (longitudes of pericenter) and n« (mean motions), where nfa^ = G(M C + rrii). The 
index i runs from the inner to the outer planet, and sets of N values can be treated as 
vectors. We assume that the parameters of the planetary system are such that the mean- 
motion resonances can be neglected. So long as the eccentricities and inclinations of the 
orbits remain small, the secular evolution of the system can be investigated using Laplace- 
Lagrange theory (Brouwer & Clemence 1961; Murray & Dermott 1999). Specifically, with 
the change of variables 

hi = Ei sin wt, ki = EiCosWi (1) 

the equations that describe the evolution of eccentricities and the motions of the pericenters 
of the orbits are 

N N 

hi = ^ ^ &jj kj ki = ^ ] A^jjhj. (2) 
j=i j=i 

The matrix A is completely determined by the masses and semimajor axes of the planets 
and the mass of the central star. In this approximation, eccentricities £j and longitudes of 
pericenter Wi are not coupled to inclinations Ii and longitudes of ascending nodes Jlj. The 
planets can be treated as coupled harmonic oscillators, and the time evolution of hi and 
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ki is a superposition of oscillations with frequencies determined by the eigenvalues of the 
matrix A and with amplitudes determined by the initial conditions and the eigenvectors of 
A. Equations similar to (2) describe evolution of inclinations and longitudes of the ascending 
nodes. 

More specifically, the matrix A is (Brouwer & Clemence 1961) 



1 m 



} N 

Am = —Hi — -^——aijaijb 3 / 2 (aij), (3) 



where 



and 



(Xij = min — , — , otij = min — , 1 (4) 



cos(ra)ofo; 



s v ; 7T 7 (1 -2a cos a; + a 2 ) 5 ' v ; 

are Laplace coefficients. With being the eigenvalues of A and e M , // = 1...N being its 
eigenvectors, the solutions of equations (2) are 

N 

h i(t) = ^e^S^smlu^t - t ) (6) 

N 

k i(t) = ^e^S^cosluj^t - t ) +/?„]; (7) 



ei (o = yfhm+m)- (8) 



The amplitude coefficients and the phases (5^ are determined by the initial conditions 
hi(t Q ) and fcj(t ). 

The initial conditions and the eigenvectors can be weighted using Wi = ai^fm^rii. 

Hi(t ) = hi(t )wi] Ki(t ) = ki(t Q )wi] E^i = e^Wi. (9) 

It can be shown that the set of vectors is orthogonal and that the weighted eccentricity 
components Hi and Ki satisfy equations similar to (2) with an orthogonal matrix A'. 

The mean square eccentricity (ef (t)) = {hf{t) + kf(t)) of each planet provides a time- 
independent measure of the level of excitation. This value can be computed from the eigen- 
vectors of the matrix A and the eccentricity components for all planets at some time t : 

m)=t (^ (H(f0> - E ;' 2+ E ( y° ),E ' i>2 ^ no) 



- 5 - 



Mean square eccentricities can be used to quantify the response of the system to an 
external perturbation - for example, another star passing not far from the planetary system. 
Before and after the perturbation, when the perturber is far enough so that its potential 
can be neglected compared to the potential of the central star and the planets, the plane- 
tary system evolves according to secular perturbation theory, as long as eccentricities and 
inclinations remain small. In this case, the mean square eccentricities after the perturbation 
are 

/J(t\\ TTf* ^2 (mo) + 5H) • E,) 2 + ((K(t ) + SK) ■ E,) 2 

\ £ iW/ a /(er~L^W (E -E) 2 ' [ f 

fx 

Here Shi an d Ski are the changes in the eccentricity components due to the perturbation 
that are weighted similarly to (9) to obtain SHi and SKi. 

Below we focus on excitation of eccentricities of planets on initially circular orbits, so 
that H(t„) = and K(t ) = 0. 



3. Excitation of eccentricities by an external perturber 

In this Section we calculate the response of the planets to an external perturber, which is 
quantified by the values Shi and Ski in the notation of the previous Section. This response can 
be obtained by numerical integration of the equations of motion in an arbitrary geometry of 
the perturbation but for simplicity we consider the two-dimensional case, when the perturber, 
the central star and the planets are coplanar at all times. We also neglect the interaction 
between the planets during the perturbation. This is justified as long as the characteristic 
timescale of the perturbation is much shorter than the timescales of secular evolution, which 
we estimate as inverse secular eigenfrequencies of the system P seC uiar ~ Porbitai,iM c /mi. 

The perturber has a mass M e and when it is far from the system it moves with velocity 
v e and impact parameter p relative to the central star (Figure 1). The equations of motion 
of the external perturber and of the planet rrii relative to the central star are 

,(M C + M e )r e _ 



—G- 



f e 



h = . G ^i-G M ^'-^-G^. (12) 
rf \Ti — r e \ 6 r s e 

Here the contribution to the potential from the planets is neglected during the perturbation; 
the planets move as test particles in the potential field of the central star and the external 
perturber. The trajectory of the perturber relative to the central star is a hyperbola. We 
restrict ourselves to the cases when the perturber passes outside the planetary system, so 
that the distance of the closest approach of the perturber is > a at. 
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The excited eccentricities remain the same if the parameters are rescaled so that a — > 
na, p — > up, v e — > K" 1//2 f e because the timescale of the perturbation r e = p/t> e scales as 
r e — > /t 3//2 r e , in the same way as does the orbital period of the planet. This can be seen from 
equations (12) which are invariant to such change of variables. 

As a result of the interaction with the central star the perturber deviates from its original 
trajectory by an angle 

V p 

y = it — 2 arctan —. - - e — — — . (13) 
A G(M C + M e ) v ; 

A perturbing star that moves with a velocity 



v. » J ^k±m (14) 

proceeds almost along a straight line with almost constant velocity. In this case (which we 
refer to as a "fast perturber") some simple approximate solutions of equations (12) can be 
obtained analytically. In Section 3.1 we present analytical estimates and numerical solutions 
for the case of the fast perturber and in Section 3.2 discuss the general case when the 
perturber does not satisfy condition (14). 



3.1. Fast perturber 

In the solar neighborhood, the velocity dispersion a ~ 40 km sec -1 , and given the 
stellar density of about n* = 0.05 pc -3 , about half of all solar-type stars have experienced 
an encounter within 500 AU at least once during the last 5 Gyr (the Sun's current age). The 
average encounter velocity is Aa/^fn in the approximation that the distribution function is 
isotropic and Maxwellian. All these encounters are fast by our definition (14), since in this 
case x ( e Q- 13) is only a few xl0~ 3 rad. 

Impulse approximation 

If the perturber is fast and if in addition the orbital period of the planet is much longer 
than the characteristic timescale of the perturbation r e = p/v e , then the displacement of the 
planet and the force from the central star can be neglected during the perturbation ( "impulse 
approximation"). Relative to the central star, the planet effectively receives a kick in velocity 
directed perpendicular to the trajectory of the perturber: 

2GM e ajcos ifi 
Aum = , 15 

V e p p-Qi COS (fi 
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where (pi is the angle between and r e at the instant of closest approach of the perturber 
to the central star. If the planet was initially on a circular orbit, then after the impact it 
proceeds on an elliptical orbit with eccentricity 



2GM e 



GM n 



V 



- x /cos 2 ^(l + 3sinV J ) + 



(16) 



For p ^> di this value reaches its maximum at (p max such that cos(2ip max ) = 1/3 and its 
minimum at (p m i n = ±7r/2. In the impulse approximation, the excited eccentricity does 
not depend on the direction of motion of the planet relative to the perturber (prograde, 
retrograde) - the effect of the direction is only important when the orbital motion of the 
planet during the perturbation is significant. 

Secular approximation 

In the opposite limit the orbital period of the planet is short compared to the charac- 
teristic timescale of the perturbation r e = p/v e ("secular approximation"). In this case the 
perturbation can be averaged over the orbital period, and thus we expect that the excited 
eccentricity is independent of the direction of motion of the planet relative to the perturber 
(prograde or retrograde) or the orbital phase of the planet at the moment of closest approach. 
The part of the disturbing function of the planet % due to the external perturber is (following 
Murray & Dermott 1999) 



R 



GM e 



i,ext 



GM e 



We expand this expression in powers of di/p: 



R 



GM e . 



i,ext 



£ 

1=2 



Pi(cosip), 



(17) 



(18) 



where ip is the angle between the vectors r» and r e . For simplicity we only consider fast 
perturbers, so that they are not significantly deflected by the interaction with the central 
star. We take the first terms (I = 2, 3) of the expansion (18) and compute the time average 
of the disturbing function over one orbit of the planet: 



{Ri,l=2) = — COS f e 

p \p 



1 ^ 1 ^ 

7 + ^ + ^>< 2 ^- 2 /e) + 0(4) 



(Ri,l=3) - 



GM e fa 



p 



— cos f e 
p 



15 

—Eicosfa - fe) +0(e1) 
16 



(19) 
(20) 



Here f e (t) describes the motion of the perturber and is a slow (compared to the orbital 
motion of the planet) function of time. For a perturber that moves almost along a straight 
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line f e is given by 



sin/ c (f) = 



v e {t - t a ) 



y/j?+V*(t-t a )*' 



(21) 



t a being the time of the closest approach. 



We substitute expressions (19)-(20) into the equations for the time evolution of eccen- 
tricities and positions of pericenters: 



ft.-+J-*W k- 

Hid* OKi 



1 d{Rj 
riiaf dhi 



(22) 



Thus, in the presence of a coplanar external perturber the eccentricity components evolve 
according to 



, GM e 3 

hi ~ — • COS J e 



p 3 rii 

h ~ COS f e 

V n i 



15 3 

/<•, ( — COs2/ e + - 



15 15 a, 2 

~7~hi sm 2j e — — COS / e 

4 16 p 



hi 



15 3\ 15 . 15 Oj . 

— cos 2/e + - + — fcj sm 2/ e - — — cos / e sm / e 
4 4/4 16 p 



(23) 
(24) 



where again we have neglected interplanetary interactions during the perturbation. These 
are approximate equations in which terms on the order of 



GM, 



p J n 



. 2 a a 
3.. ■ \ £ >^> 6 - 



p2 p 



(25) 



are neglected. 



For a planet on an initially circular orbit, in equations (23)-(24) the terms proportional 
to k and h can be neglected as long as k, h a/p during the event. In this case equations 
(23)-(24) can be integrated to obtain 



5 hi 



5GM e 



di at 



4 pv e V GM C p 2 



, Ski — 0, Sei = \5hi\. 



(26) 



This result can also be obtained from general formulas presented by Heggie & Rasio (1996) 
and Kobayashi & Ida (2001) 2 in the limit when the perturber is coplanar and moves very 
fast. 



2 Due to a typographic error there is a factor of e» 3 missing in equations (46) and (47) of Kobayashi & 
Ida (2001). 
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Intermediate regime 

If a planetary system has a large range of semimajor axes, it is plausible that the same 
passing star acts as an "impulse" perturber for the outer planets and as a "secular" perturber 
for the inner planets. Below we consider a case of a perturber with M e = M Q that has impact 
parameter p = 500 AU and velocity v e = 80 km sec -1 . In this case the deflection angle of 
the perturber x = 1 x 10~ 3 . We numerically integrated the equations of motion (12) for 
M c = M Q and obtained values of the excited eccentricities for a large range of semimajor 
axes (from 1 AU to 300 AU) and for several values of the orbital phase at closest approach 
if (cf. Kobayashi & Ida 2001). The results are presented in Figure 2 (top) together with the 
analytical approximations (26) and (16). 

The secular approximation (26) is shown with a dashed line. It is applicable for planets 
with orbital periods P or utai ^ T e, or in terms of semimajor axes a C 7 AU for the encounter 
with p = 500 AU and v e = 80 km sec -1 , and the analytical formula (26) describes the 
excitation very well throughout its region of applicability. In the intermediate regime (a of 
about tens of AU) the excited eccentricity depends both on the orbital phase at the closest 
approach of the perturber and on the direction of the motion of the planet relative to the 
perturber, and the average excited eccentricity is substantially larger for prograde planets. 
For the slowest planets (large a) the impulse approximation (16) works reasonably well and 
the excited eccentricity only depends on the phase, not on the direction, in agreement with 
equation (16) and related discussion. The discrepancy between the theoretical prediction 
and the numerical calculation is due to neglecting terms of the next order in a/p in (16) 
which become significant at a ~ 100 AU. 



Many or most stars are born in clusters. Open clusters are relatively short-lived (~ 10 8 
years) systems with velocity dispersions of 1 — 3 km sec -1 that reflect their small mass and 
with central densities about 



where r c is the cluster core radius (King 1962; Lynga 1982). In such environments about 
half of the stars experience an encounter with impact parameter p < 400 AU once during 
the lifetime of the cluster, but condition (14) no longer holds, and deflection of the perturber 
becomes important. 

In Figure 2 (bottom) we show the excited eccentricities calculated numerically from 



3.2. Slow perturber 




(27) 
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equations (12) as a function of semimajor axis in an impact with p = 400 AU and v e — 5 
km sec -1 . In this case the periastron distance of the perturber is 335 AU. Our secular 
approximation (26) is no longer valid because the perturber does not move along a straight 
line anymore, but generalized analytical solutions in the secular approximation were obtained 
by Heggie & Rasio (1996) and Kobayashi & Ida (2001) 3 : 



_ 15 M e /M c fdi\ 5/2 1 



8 y/l + M e /M c V P J e H (e 2 H - I) 5 / 4 
where 



A arccos ( -— ) + ^-iL (1 + 2e' H) 



(28) 



is the eccentricity of the hyperbolic motion of the perturber. Approximation (28) is plotted 
with a dashed line in Figure 2 (bottom). 

Typical eccentricities excited during a slow interaction are one or two orders of magni- 
tude higher than those excited by a fast perturber passing at similar distances. The expected 
number of encounters within a given impact parameter is similar in the solar neighborhood 
and in a rich open cluster; the large density of stars in the latter compensates for the short 
life-time and slow velocities. The result is that eccentricities acquired at the birthplace of 
stars during the first 10 8 years are much higher than those that can be generated during the 
subsequent evolution in a low-density environment. In what follows we will thus concentrate 
on slow encounters. 



4. Propagation of the disturbance via secular interactions 

After the perturber is gone, the excited eccentricities are redistributed in the system due 
to the secular interactions described in Section 2. As a first example, we consider a planetary 
system around a solar mass star consisting of 10 planets with their semimajor axes following 
a geometric progression from 1 AU to 100 AU, with / = Oj+i/aj = 1.67. The masses increase 
as a geometric progression, too, with the innermost planet having the mass of 0.1 Mj up i ter 
and with g = mi+i/rrii = 1.67. The outer planet has then a mass of itln = 10Mj np j ter . 
The planets were on circular orbits before the perturbation, and the orbital phases at the 
moment of the closest approach (p,i were selected randomly. The perturber is in prograde 
motion relative to the planets. 



3 Duc to a typographic error there is a factor of e„ 3 missing in equations (46) and (47) of Kobayashi & 
Ida (2001). 



- 11 - 



The evolution of the eccentricity of the innermost planet is shown in Figure 3 (top). If 
there were no other planets in the system, a perturber passing with an impact parameter of 
400 AU and with velocity v e = 5 km sec -1 would only excite an eccentricity of 1.75 x 10~ 7 . 
However, secular interactions with other planets bring the eccentricity up to a maximum of 
0.015. The rms eccentricity can be found from equation (11) to be \J (£?(£)) — 0.008. 

In a more extreme example given in Figure 3, bottom, the spacing between the planets 
is the same but the masses increase from mi = M Earth to itln = Mj up i ter with g = 1.90. In 
this case secular interactions in the extended planetary system bring the eccentricity of the 
innermost planet up to 0.25, with an rms value of yj (s\{t)) — 0.13. 

In both examples, the presence of outer planets made the perturbation much more 
efficient in exciting the eccentricity of the inner planet. We define for the inner planet an ex- 
citation efficiency parameter which is the ratio of the rms eccentricity after the perturbation 
to the value of eccentricity that would have been excited if there were no other planets: 

X = \J {e\ (after perturbation) )/e 1 (excited if no interactions). (30) 

In the cases considered above, the excitation efficiency for the inner planet is X ~ 4.6 x 10 4 
(Figure 3 top) and X ~ 7.2 x 10 5 (Figure 3 bottom). 

Another important parameter is the ratio of the eccentricities of the innermost and the 
outermost planets. If this ratio is <SC 1, the eccentricity of the innermost planet cannot be 
excited significantly without ejecting the outer planet. We define a propagation efficiency 
parameter: 

P = \J {e\ (after perturbation))/ \J (e 2 N (after perturbation)). (31) 

In the case of our example system in Figure 3 (top), the rms eccentricity of the outermost 
planet y/(e 2 N (t)) ~ 0.029, so the propagation efficiency is P ~ 0.28. In the second example 
(Figure 3 bottom), y/(e 2 N (t)) ~ 0.029 and P ~ 4.4. So in fact, the final rms eccentricity of 
the innermost planet is the largest in this system even though this planet is least affected 
by the direct interaction with the perturber. 

Both efficiency parameters are quite high for our example systems, but they depend 
significantly on the mass distribution. The latter can be parametrized by the surface density 
that is obtained by smearing out the mass of all planets in such a way that the mass in the 
annulus between a« and Oj+i is rrii. The resulting smeared-out surface density is a power 
law as a function of the distance from the central star, and the power index is related to the 
ratios of semimajor axes and masses: 



E(r) oc r~ g , q = 2- ln(g) / ln(/). 



(32) 
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We now proceed to investigate how the efficiencies X and P depend on the parameters 
of the planetary system. We calculated rms eccentricities for 16,000 simulated systems, 
with the total number of planets N ranging from 5 to 20. The planets follow a geometric 
progression in semimajor axes and masses quantified by the parameters / = Oj+i/aj and 
g = TOj+i/mj. All planets move in the same direction and are on initially circular orbits, 
and the perturber always has the same parameters (M e = M & , p = 400 AU, v e — 5 km 
sec -1 ). For each N, the parameter / was selected randomly, with the maximum size of the 
system a N restricted to be in the range 50-200 AU. For each pair N, f the parameter g was 
selected randomly, but the ratio of the maximum mass to the minimum mass was not to 
exceed 5M Jupiter /M Earth ~ 1600. For each set N, f,g we selected random sets of phases ^ 
and calculated the resulting rms eccentricities for both prograde and retrograde motions of 
the perturber relative to the planets. The innermost planet was placed at 1 AU from the 
central star in all systems, so that in all systems 



In Figure 4 we plot the rms eccentricity of the innermost planet and the propagation 
efficiency for this planet as a function of the surface density index q. Most notably, both 
the excited eccentricity and the propagation efficiency are high in systems with flat mass 
distributions q < 1. Prograde encounters are, on average, about ten times more efficient in 
exciting the eccentricity than retrograde encounters. It can be seen from Figure 4 that high 
eccentricities can be excited for flat density distributions. The Laplace-Lagrange secular 
perturbation theory cannot be used for eccentricities approaching unity, and therefore we do 
not show systems with \J (ef (£)) > 1 (numerical integrations showing the effects of secular 
interactions of high-eccentricity planets are discussed by Terquem & Papaloizou 2002). We 
also do not show systems in which the outer planet is ejected as a result of the direct 
interaction with the perturber (this occurred in about 9% of systems with sizes in excess of 
100 AU; all planets with semimajor axes less than 100 AU remained bound to the central 
star). 

Our simulated systems occupy only a narrow locus on the q vs efficiency diagrams, but 
this effect is artificial and is introduced by the restrictions we imposed on the simulated 
systems. In particular, we allowed only a narrow range of total sizes of the systems (50-200 
AU). It is natural to assume that the excited eccentricities also depend on the total size of 
the system and, regardless of the mass distribution, high eccentricity cannot be excited in 
small systems. We therefore added to our simulation 11,000 more systems with sizes in the 
range 10 AU < a N < 50 AU and studied the dependence of the excited eccentricity on the 
total size of the system. For the encounter that we studied (p = 400 AU, v e — 5 km sec -1 ), 
the excited eccentricity is almost insensitive to size for 50 AU < a N < 200 AU, but drops 




(33) 
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significantly for smaller systems (Figure 5). 

Within the Laplace-Lagrange theory the final eccentricities do not depend on the overall 
mass scaling, but only on the ratios of the masses of the planets. The stability of the system, 
however, depends on the actual masses. Gladman (1993) gives a stability criterion for two- 
planet systems; if this criterion is applied to our simulations, the planetary system is stable 
if / — 1 > 3(m max /M c ) 1 / 3 where m max is the maximum mass of the planets in the system and 
M c is the mass of the host star. All simulations with 5 < iV < 20 and the total sizes 50 — 200 
AU are stable by this criterion if m max < OAMj. If higher masses are allowed, the stability 
criterion limits the number of planets in the system. For example, if m ma x = 5Mj, then 
only systems with N < 13 can be stable. The actual requirements for the long-term stability 
of multi-planet systems are likely to be somewhat more severe than Gladman's criterion 
suggests. Both excitation and propagation efficiencies are insensitive to the number of the 
planets in the system, and excluding unstable systems does not affect any results of our 
analysis. 

We also performed the same simulations for fast encounters (p = 400 AU, v e = 80 
km sec^ 1 ). The excited eccentricities are typically about two magnitudes lower than in the 
case of slow encounters, in agreement with Figure 3, but the propagation efficiencies are 
very similar to those shown in Figure 4. The main factors governing the efficiency of the 
eccentricity excitation and propagation are the total size of the planetary system and its 
large-scale surface density distribution, with efficient eccentricity excitation and propagation 
in systems with q < 1. This conclusion is valid both for fast and slow encounters that we 
examined. 



5. Discussion 

We showed that encounters with nearby stars can excite the eccentricities of the outer 
planets, which can then propagate to the inner planets. Eccentricity propagation is efficient 
for flat mass distributions of the planetary system. We now discuss the typical values of 
eccentricities excited by this mechanism in different environments. 

In the solar neighborhood, all encounters are fast by criterion (14), and the impulse 
approximation works fairly well for planets at radii of a few tens of AU (Figure 2 top). After 
the eccentricity excitation propagates to the innermost planet, its rms eccentricity becomes 

rms( £l ) = 0.0025P ( - T ) V ^ Y'* ( P x Y* = 0.01 -0.1 (34) 

K J V80kmsec-V vl00 AU/ V500 AU/ V ; 

for P = 4 — 40. We have averaged over all phases of closest approach in equation (16). 
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In an open cluster, the secular approximation (28) sets a lower bound on the excited 
eccentricity (Figure 2 bottom), and for planets at radii larger than 50 AU the average excited 
eccentricities are about a factor of 3 larger than those given by the secular approximation 
(28). The rms eccentricity of the inner planet is then 

rms( £l ) = 0.053P ( V * Y f^-f 2 (— > 0.2 (35) 
K J VSkmsec -1 / V100AU/ V400 AU/ v ; 

for P > 4. 

Typical eccentricities excited in dense open clusters are a factor of 20 greater than those 
that can be excited by field stars in the solar neighborhood. In the simulations shown in 
Figure 4 propagation efficiencies of P = 1 — 100 were achieved for flat mass distributions 
(q < 1). If such propagation efficiencies are common, the high eccentricities typical of those 
observed in extrasolar planets can be produced by close encounters in birth clusters during 
the ~ 10 8 years after the stars and planetary systems are formed. Future observations of 
systems containing extrasolar planets may determine whether there is enough mass in the 
outer parts of these systems for the suggested mechanism to be efficient. The observation 
times required to detect radial velocity signatures from companions at radii > 50 AU are 
prohibitively long, but such objects can be detected directly by infrared imaging using the 
Spitzer telescope or by high-contrast imaging. The suggested mechanism also operates if the 
mass in the outer parts of the system is in the form of planetesimals. The detectability of 
such disks depend on their physical conditions. 

The probability that the closest encounter during the life time of the system r is p or 
smaller is given by Poisson's formula: 

n(< P ) = i_ e -™.w (36) 

Here n* is the number density of the stars that act as perturbers. Therefore, about 10% 
of all stars experience very close encounters (within 200 AU). In these cases the excited 
eccentricity can be very high even for field encounters in the solar neighborhood and is very 
sensitive to the size of the system and the phase of the outer planet at closest approach of 
the perturber. 

6. Conclusions 

The eccentricities of most extrasolar planets are much larger than those of giant planets 
in the Solar System. This observation requires explanation because planets that form from 
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a protoplanetary disk are expected to have nearly circular orbits. In this paper we explore 
the possibility that close stellar encounters can excite planetary eccentricities. 

Encounters with nearby stars can be dynamically important for the outer parts of cir- 
cumstellar disks and planetary systems. Kobayashi & Ida (2001) studied pumping up of 
eccentricities and inclinations of planetesimals in a close stellar encounter. Kalas et al. 
(2001) and Larwood & Kalas (2001) suggested that encounters with nearby stars can be 
responsible for asymmetries in the j3 Pictoris disk. Hurley & Shara (2002) found that most 
planets are stripped from their host stars by encounters with neighbor stars in the dense 
centers of globular clusters. 

In the solar neighborhood or an in open cluster the eccentricities of short-period plan- 
ets cannot be excited directly by perturbations from nearby stars to any significant values 
because such planets are tightly bound to their host star. We suggest here that if there are 
other planets or a significant mass in smaller bodies in the same system on larger orbits, the 
outer planets can be excited to high eccentricities and then transfer this excitation to inner 
planets through secular perturbations. 

The extent to which the eccentricity of the innermost planets can be excited depends 
on the total size of the planetary system and on its smeared-out surface density distribution. 
If the smeared-out surface density distribution is parametrized by a power law, S(r) oc r~ q , 
then the proposed mechanism is important for systems with q < 1 and with total size 
^at > 50 AU. At least some gas disks around young stars show flat density distributions 
with q ~ 1 (Dutrey et al. 1996; Wilner et al. 2000) and extend out to tens or hundreds of 
AU. Theoretical models of steady-state protoplanetary disks also suggest flat surface density 
distributions (Bell et al. 1997). It appears that the maximum lifetime of such gas disks is 
10 6 — 10 7 years (Haisch et al. 2001); some of this material may dissipate, but much or most 
of it may survive as planets or large planetesimals. Outward migration of planets can also 
yield massive planets at a few tens of AU (Veras & Armitage 2004). It has been suggested 
that the asymmetries in the dust debris disks in some systems may be due to unseen planets 
at radii 30 - 60 AU (Ozernoy et al. 2000; Wilner et al. 2002). 

If the mass distribution is flat and the planetary system is large, the inner planets can 
be excited to high eccentricities (0.1 — 1, eq. 35) due to encounters in rich birth clusters. 
The typical eccentricities that can be excited by the described mechanism from field star 
encounters in the solar neighborhood are smaller (0.01—0.1, eq. 34). A small fraction 
of stars (10%) experience very close encounters with field stars (within 200 AU) in which 
excited eccentricities can be much higher. 

All the calculations were done under a simplifying assumption that the perturber is 
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moving in the same plane as the planets. Kobayashi & Ida (2001) argued that the eccen- 
tricities excited by a passing star are insensitive to the inclination of the perturber relative 
to the planetary disk and lie in the range determined by coplanar prograde and retrograde 
encounters. Therefore, our conclusions about the eccentricity excitation and propagation 
are likely to be independent of the inclination of the perturber. In the case when the trajec- 
tory of the perturber is inclined relative to the plane of the planetary system, inclinations 
of the outer planets in the system can be excited and then propagate inward, similarly to 
propagation of eccentricities that we discussed in this paper. Warps in circumstellar disks 
like those seen in (3 Pictoris (Kalas & Jewitt 1995; Wahhaj et al. 2003) may be signatures 
of inclination waves propagating through the system after the initial perturbation. 

In this work our focus was to find the conditions for efficient propagation of eccentricity 
disturbances. In our simulations, the masses and semimajor axes of successive planets follow 
geometric progressions. Clearly, these special configurations do not represent the full range 
of properties of real systems. Therefore, we do not discuss the distribution of eccentricities 
resulting from our simulations in relation to the observed distribution. 

We would like to thank Marc Kuchner and Roman Rafikov for useful discussions and 
the referee for helpful comments on the manuscript. ST acknowledges the support of NASA 
grants NAG5-10456 and NNG04GH44G. 
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Fig. 1. — Geometry of the encounter in the reference frame of the host star M c : the planetary 
system orbits the star M c and the passing star M e is a perturber moving on a hyperbolic 
trajectory. v e is the velocity of the perturber far from M c , and p is its impact parameter. 
In this Figure, the motion of the planet relative to the perturber is prograde if it is moving 
counter-clockwise and retrograde otherwise. 
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Fig. 2. — Top: Excited eccentricity as a function of the semimajor axis a of the planet for a 
fast encounter (with parameters typical for the solar neighborhood). Solid lines are results 
of the direct integration of the equations of motion for two values of the orbital phase tp at 
the moment of closest approach of the perturber, and for two directions of the motion of 
the planet relative to the motion of the perturber (prograde and retrograde). Two extreme 
values of the phase (defined after eq. 16) are chosen, so that the results for other phases lie 
in between the uppermost and the lowermost curves. The dashed lines show the analytical 
approximations in the secular regime (eq. 26). The dotted line is the impulse approximation 
(16) plotted for <p = ip ma x] it represents the theoretical maximum of the excited eccentricity. 
The results of the numerical integration for ip = (p max deviate from the dotted line because 
the assumption a p breaks down for large a. 

Bottom: Same as above but for a slow encounter with parameters typical of rich open clusters. 
The dashed line is the secular approximation by Heggie & Rasio (1996) and Kobayashi & 
Ida (2001). 
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Fig. 3. — Evolution of the eccentricity of the innermost planet due to secular interactions 
with other planets in two different systems. Top: masses of the planets range from 0.1 
Mj U pi ter to 10 Mj U piter- Bottom: masses of the planets range from MEarth to Mj upiter . The 
parameters of the systems and of the perturbation are specified in the top left corner of 
each diagram. If there were no outer planets the eccentricity of the innermost planet would 
be just 1.75 x 10~ 7 as a result of the direct interaction with the perturber. In both cases 
the secular evolution is dominated by perturbations from the two outermost planets; the 
characteristic timescales are P seC uiar,i ~ PorUtai,iM c / rrii with i = N producing the slow large 
amplitude variation and i = N — 1 producing smaller amplitude variation. 
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Fig. 4. — Excitation and propagation efficiencies for a slow perturbation as a function of 
the smeared-out surface density distribution parameter. Prograde encounters are shown in 
black, and retrograde ones are shown in red; the parameters of the simulation are listed in 
the top right corner of both panels. In our ensemble of simulated systems all innermost 
planets are located at 1 AU and the parameters of the perturber are the same; therefore, the 
excitation efficiency is directly proportional to the rms eccentricity of the inner planet: X = 
a/ (ef (£))/1.75 x 1CT 7 . Systems in which the outer planet is stripped by the perturbation are 
not shown, and systems with very high calculated eccentricities >1 are shown as \J (£i(t)} = 
1. 
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Fig. 5. — Dependence of excitation and propagation efficiencies on the total size of the 
system. Only prograde encounters are shown, and we zoomed on the high-efficiency part 
of Figure 4. All systems are divided in four groups according to their size, and the median 
efficiency parameters are shown as a function of the mass distribution parameter q for each 
of the three size groups. The total number of planets in this simulation varied between 5 and 
20, the innermost planet was always at ai = 1 AU and the parameters of the impact were the 
same: p = 400 AU, v e = 5 km sec -1 , M c = M e = M & . The excited eccentricity is insensitive 
to the total size of the system as long as 50 < a N < 200, but drops significantly in 
systems smaller than 50 AU. 



